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Abstract -We investigate the effect of mobility on the response of coupled oscillators to a sub- 
threshold external signal in metapopulation networks, wherein each node represents a subpop- 
ulation with overdamped bistable oscillators that can randomly diffuse between nodes. With 
increasing mobility rate, the oscillators undergo transitions from intrawell to interwell motion, 
demonstrating clearly mobility-enhanced signal amplification. Moreover, the response shows non- 
monotonic dependence on the mobility rate, i.e., a maximal gain occurs at a moderate level of 
mobility. This interesting phenomenon is robust against variations in the overall density, net- 
work size, as well as network topology. In addition, a simple mean-field analysis is carried out to 
qualitatively illustrate the simulation results. 
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Introduction. — Enhancing the collective response to 
a weak input signal is an important and challenging prob- 
lem in a variety of fields, not only in that of traditional 
signal processing but also in those such as particle physics 
[1] , gravitational wave search [2] , and medical science [3ll4] . 
Many natural and artificial information-processing sys- 
tems are connected together to form functional networks 
and spontaneously adjust their internal machinery to en- 
hance the sensitivity to external signals. For instance, 
cells and microorganisms respond to changes in external 
environment by means of an interconnected network of 
receptors, messengers, protein kinases and other signal- 
ing molecules [5H7]- One of the most intriguing part of 
these phenomena was amplified signal response. It has 
been shown that random fluctuations can enhance the re- 
sponse to weak periodic driving, as observed in many dif- 
ferent physical, chemical and biological scenarios [5HTU]. 
Recently, amplified signal response in complex network 
of coupled oscillators has drawn considerable attention 
[71 lllftT5] . It has been found that the weak external sig- 
nals can be amplified by the heterogeneity in degree of 
the network [TlHHIIl] , adaptive coupling weights between 
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the signal node and its neighbors [13| . neuronal diversity 
on complex networks |14j and a multilayer feedforward 
network |15| . A one-body theory which gave analytic ex- 
pressions of the gain and the degree of the unit with the 
maximum response to the input signal in terms of the cou- 
pling strength was also developed [TB) . 

Nevertheless, previous studies on signal response in 
complex networks only deal with the case of immobile ele- 
ments and each network node is occupied by one single el- 
ement. Very recently, the metapopulation network model 
|17] . which incorporates subpopulation in the node, mo- 
bility over the nodes, and a complex network topology, 
has attracted much attention. This model has been suc- 
cessfully exploited in different contexts, such as epidemic 
spreading P"5H2D] . biological pattern formation pTlf^ . 
chemical reactions [53], population evolution and 
many other spatially distributed systems [ISJI^S]. It is 
shown that the mobility and the density of the individu- 
als could have drastic impact on the emergence of collec- 
tive behaviors in general [T7l[27], and particularly, mobil- 
ity induced and tuned synchronization of coupled oscilla- 
tors have been reported pSMO] . Therefore, one may ask: 
How would the mobility influence the signal response in 
metapopulation networks of coupled oscillators? 
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In the present work, we consider a metapopulation net- 
worked system wherein each node is occupied by any num- 
ber of overdamped bistable oscillators, subject to a sub- 
threshold external signal. By stochastic simulations of the 
involved dynamical reaction-diffusion processes, we find 
that, the signal response exhibits a clear-cut maximum at 
an optimal level of mobility rate. Furthermore, we show 
this nontrivial phenomenon is robust to the density and 
network size as well as different network topologies. A 
simple mean-field (MF) analysis is given to help us under- 
stand the simulation results. 

Model description. — Wc consider a system of M 
individuals distributed in distinct subpopulations la- 
beled /i, each corresponding to a network node, and as- 
sume that the number of individuals in node ^ is N^j^, 
satisfying X]^=i -^a' ~ ■ Thus, the density p of the 
metapopulation is given by p = M/N. Individuals inside 
each subpopulation run through the paradigmatic bistable 
oscillators, and the dynamics of the ith-oscillator located 
in the /ith node is described by: 

- X,) + Asin{ujt) (1) 



Here Xi (i = 1, M) is the state variable of the i-th unit 
at time t, and C is the coupling strength. A and uj are the 
amplitude and frequency of the external periodic forcing, 
respectively. 

The above equation actually defines the "reaction" pro- 
cess that governs the overdamped motion of a Brownian 
particle in a double- well potential with subthreshold peri- 
odic forcing. We now assume that the particles can diffuse 
randomly among the nodes. The system evolves in time 
according to the following rules [IT]. We introduce a dis- 
crete time step r representing the fixed time scale of the 
process. The reaction and diffusion rates are therefore 
converted into probabilities. In the reaction step, all the 
particles are updated in parallel according to EqUJ Af- 
ter that, diffusions take place by allowing each particle to 
move into a randomly chosen neighboring node with prob- 
ability Vt, where V denotes the mobility rate. If not oth- 
erwise specified, the parameters are A^ — 1000, r = 0.001 
and C ~ 0.005. We choose the mobility rate V as the 
main control parameter. Each simulation plot is obtained 
via averaging over 50 independent runs. 

Simulation results. — To begin, we consider scale- 
free networks generated by using the Barabasi-Albcrt 
(BA) model [31] with power-law degree distribution 
p{k) k^^ and average degree (k) = 6. We fix p = 10 
and vary V to investigate how the state variables of oscilla- 
tors evolve in time. Initially, the oscillators are randomly 
distributed among the nodes and homogeneously assigned 
the initial position at ±1, which are the two minima of 
the double-well potential. The signal is considered to be 
subthreshold, i.e., it does not suffice to induce the oscil- 
lators to jump between the two minima in the absence of 



diffusion. For small V, the oscillators are separated into 
two distinct subsets: Some oscillate around the minimum 
1 and the others around -1, depending on their initial con- 
ditions. As examples, typical time series of Xi(t) for sev- 
eral randomly-chosen oscillators arc plotted in Fig. 1(1) for 
V — 2.5x10"''. For moderate mobility rate F = 2.5x10"^ 
as shown in Fig. 1(11), jumps between the two stable wells 
tend to occur periodically in time, driven by the periodic 
force. However, for large V, say V = 0.25 in Fig.l(III), the 
oscillators are synchronized but all confined into one single 
well. Therefore, we observe interesting mobility-induced 
transitions from intrawell to interwell motion, and then to 
synchronized single-well motion. We note that the final 
state for large V may depend on the amplitude or the fre- 
quency of the external signal. For instance, for A ~ 0.38 
(w = 0.015), the transition from intrawell to internal well 
motion still occurs when V increases from small values, 
but the oscillators finally oscillate separately around the 
two mimima if V is large enough(not shown). Wc also 
note that this nontrivial transition and rc-entrancc phe- 
nomenon is observed if one fix V and let p change. In 
this latter case, the oscillators will remain in the well they 
initially are no matter how large p is. 
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Fig. 1: (Color online) Typical time evolutions of the state 
variables Xi{t) for several randomly-chosen oscillators sd, V — 
2.5 X 10"^ (upper panel), V = 2.5 x lO"'* (middle panel) and 
V = 0.25 (lower panel). Time in the simulations is measured in 
units of Monte Carlo steps (MCS), where one MCS is defined 
as M reaction and diffusion attempts, M being the number of 
the total oscillators in the metapopulation network. It is shown 
that the time series undergo mobility-induced transitions from 
intrawell to interwell motion, and then to synchronized single- 
well motion. Other parameters are A'^ — 200, (k) — 6, p = 10, 
A = 0.39 and uj = 0.015. 
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To quantitatively measure the signal response, we intro- 
duce a gain factor G which is defined as follows [71[T31[T1] , 



G 



1 

AM ^*=i 



maxxi — minXi)/2 



(2) 



where [•] stands for averaging over 50 different network re- 
alizations for each value of V . Clearly, a large G means a 
larger signal response. Figsl2]plot the G as a function of 
V for fixed p = 10 with varying signal amplitude A or fre- 
quency ui. Clearly, G shows a nonmonotonic dependence 
on the mobility rate as expected from Fig.l. Note that 
for large V ^ the final stationary values of G may have two 
different values depending on uj and A in a somewhat com- 
plicated way: One approaches 1.5 for small uj and large A, 
as depicted by red solid circles and dark yellow dotted cir- 
cles in FiglUJa) and (b), and the other approaches 0.75 
for small w with small A, or for relatively large w, as also 
shown in Fig. 2. The time series shown in Figlljcorrespond 
to different mobility rates indicated by the three arrows in 
Fig. 2(b): The transitions from intrawell to interwell and 
then to synchronized single- well motion are demonstrated. 

So far, the results are all for scale- free coupled networks. 
One may wonder whether the above interesting findings 
are sensitive to the network topology or not. Thus, we 
have also performed similar studies on other types of net- 
works, such as small world networks and random networks. 
Fig El plots the dependences of G on F for a typical small- 
world network and a random network, shown by triangles 
and squares respectively. Apparently, the qualitative be- 
haviors are the same as those observed for scale-free net- 
works. The only difference is that the optimal mobility 
rate for the maximal signal response is different. 

Mean field analysis. — To get more insight into the 
aforementioned results, here we present a simple mean- 
field analysis by considering a model system consisting of 
only two nodes. This makes the problem mathematically 
tractable, while still capable of capturing the main trait 
of mobility. Assuming that the node 1 and 2 hold iVi and 
N2 oscillators respectively, we can introduce the average 
state variables Xi = ]^ X^iei and X2 = ]^I]jg2^J ' 
whose dynamics is governed by the following equations 
via averaging Eq. (1) over each subpopulation, 

Xi=Xi-^"^x^^ + A sm{ujt) + V{X2-Xi) (3a) 
X2 = X2- + A sm{ujt) + V{Xi - X2) (3b) 

The first three terms in the right-hand side of Eq. (3a) 
account for averages with respect to the counterparts in 
the Eq.(l), and the last term represents the diffusion of 
elements into and out of the node 1. Note here that 
we have chosen the characteristic time r as the discrete 
time step to integrate Eq.(l). Equation (3b) can be in- 
terpreted in a similar manner. Following the scheme 
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Fig. 2: (Color online) The gain factor G as a function of mo- 
bility rate V on BA scale-free networks at different uj for fixed 
A — 0.39 (a), (b) and (c) corresponding to different A at fixed 
OJ — 0.015 and uj — 0.02 respectively. The representative time 
evolutions at different mobility rates indicated by arrows in 
panel (b) are shown in Fig[T] 



used in Ref. |32j . one may introduce the variances of 
state variables within the node 1 and node 2, denoted by 

o"! = I^iGi - ^1)^ a^'i = 77; Eje2 - ^2)^ 
respectively. Assuming these variances to satisfy Gaus- 
sian distribution, Eqs. (3) become 

A1.2 = Ai.2(l - 3^1.2) ~Xl2 + A sin(wt) + V{X2,i - Xi,2) 

(4) 

However, this equation is yet not solvable analytically 
since we do not know the exact expressions of i7i_2- To 
proceed, we numerically calculate ai^2 by using an ensem- 
ble average along with the numerical integration of Eq.(4). 



p-3 



C. Shen et al. 




1E-5 1E-4 1E-3 0.01 0.1 1 10 100 
V 



Fig. 3: (Color online) The gain factor G as a function of 
mobility rate V in the metapopulation model, triangles and 
squares correspond to small-world network and random net- 
work respectively, both on a synthesized 1000-node network 
with (fc) = 6. 

The gain factor associated with this simple model is then 
obtained as G = 2^ Y^i=i (niaxX^ — xmnXi)/2. In Fig|31 
we give the results of G as a function of the mobility rate 
V for the two-node network, where the lines denote the re- 
sults obtained from Eq.(4) and the symbols from Eq. (1). 
Clearly, the mean field equation (4) can reproduce quali- 
tatively well the main character: There exists an optimal 
mobility rate where the gain reaches the maximum. 

Discussion and conclusion. — It is now well known 
that mobility may play constructive roles in many coupled 
systems. For example, mobility can lead to the optimal 
synchronization in two-dimensional coupled phase oscilla- 
tor model [33] and partially occupied networks [34], tune 
synchronization of integrate-and-fire oscillators |30| , influ- 
ence the synchronization pathway in metapopulations of 
mobile agents [H] and affect the epidemic threshold in 
metapopulation networks |35| , etc. Our findings here show 
that natural systems might profit from their mobility in 
order to optimize the response to an external stimulus. 
Recently, we have also found that mobility can consider- 
ably induce metapopulation coupled oscillators undergo- 
ing phase transitions from incoherent to amplitude death, 
and then to synchronized state [2S|. Such constructive 
effects of mobility in complex systems may deserve more 
and more attention in future works. 

In summary, we have studied the signal response of 
coupled bistable oscillators in metapopulation networks, 
wherein different subpopulations are connected by fluxes 
of individuals. By extensive numerical simulations, we 
show that the mobility plays nontrivial roles on the col- 
lective response of the system, by demonstrating an in- 
teresting type of mobility-enhanced signal response to an 
external periodic forcing, which is robust to the density, 
network size as well as network topology. We have also 
performed a simple mean field analysis which can qual- 
itatively reproduce the simulation results. Since many 
real- world networks, such as cellular networks, protein 
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Fig. 4; (Color online) The gain factor G as a function of mobil- 
ity rate V on two-node networks. The lines denote theoretical 
results and the symbols denote simulation ones. Other param- 
eters are A = 0.388 and p = 200. 

networks, gene networks, and so on, inevitably involve 
variances in the mobility, and their collective dynamics 
could be modeled by coupled oscillators, these results may 
find many applications in several fields of physics, neuro- 
science, and biology. Our study may provide valuable in- 
sights into the mobility-induced collective response to the 
external signals that take place in other metapopulation 
networked systems. 
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